Surface plasmons at composite surfaces with diffusive charges 
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Metal surfaces with disorder or with nanostructure modifications are studied, allowing for a lo- 
calized charge layer (CL) in addition to continuous charges (CC) in the bulk, both charges having 
a compressional or diffusive non-local response. The notorious problem of "additional boundary 
conditions" is resolved with the help of a Boltzmann equation that involves the scattering between 
the two charge types. Depending on the strength of this scattering, the oscillating charges can be 
dominantly CC or CL; the surface plasmon (SP) resonance acquires then a relatively small linewidth, 
in agreement with a large set of data. With a few parameters our model describes a large variety of 
SP dispersions corresponding to observed data. 
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Introduction. Collective electronic excitations in 
metal surfaces covered with adsorbates or nanostructures 
are of significant recent interest. Surface plasmons (SP) 
are an efficient tool for characterizing such surfaces^'^ and 
can be used as sensitive chemical sensors and biosensors^. 

A considerable amount of data on the dispersion of SPs 
has been accumulated^"— on a variety of metal surfaces, 
clean, sputtered, or covered with thin films. These com- 
posite surfaces indicate the necessity of distinguishing be- 
tween two types of charge carriers: continuous charges 
(CC) extending throughout the bulk and a charge layer 
(CL) of carriers localized at the surface. In fact, photo- 
emission data on some of these surfaces reveals the exis- 
tence of quantum well states at the surface^^. The two 
charge types are relevant also for pure metals: in alkali 
metals the charge extends beyond the neutralizing ions, 
forming a distinct layeri^. In Ag a two-component s-d 
electron system with different surface and bulk charges 
has been put forward to explain the SP dispersioni^. In 
some cases a distinct surface band is formed^, e.g. as in 
Be(OOOl), showing an acoustic plasmoniS. Further moti- 
vation for a two-type charge model comes from studies of 
the anomalous heating of cold ions observed in miniatur- 
ized Paul traps, that invoke surface charge fluctuations 
on the metallic electTodes^"—. 

Most information about the SP dispersion uj{k) is avail- 
able in the non-retarded range ujp/c 'ti k kp where k 
is the momentum parallel to the surface, ujp is the bulk 
plasma frequency, c the speed of light, and kp the Fermi 
momentum. In this range, the dispersion is parameter- 
ized a.s uj{k) = A + Bk -\- Ck^ and the limiting value 
A — Wp/v^ is well known (assuming unit background 
permittivity)^^. Considerable insight is gained by Feibel- 
man's sum rule^ relating the slope B to the centroid of 
the oscillating charge density profile Sn{z): 

g ^ Jdzz6n{z) 
2V2 Jdz6n{z) 

Here the medium is located in the half-space z < 0, and 
the onset of the dielectric function due to bulk ions is at 
z = 0. Hence the SP dispersion is negative {B < 0) if 



the oscillating charge is dominantly located outside the 
metal (z > 0) as in alkali metals^*^, while -B > if the 
oscillating charge is dominantly inside. In some cases the 
coefficient B is small and the quadratic term Ck^ domi- 
nates the dispersion^iiiii^, even at the lowest measurable 
k. From Eq.([T]), i3 = is a strong indication for a charge 
excitation that is highly localized at the surface, z = 0. 

In the present work we consider a non-local electrody- 
namic model including both CC and CL, which for small 
k involves diffusive (or compressional) terms, similar to 
hydrodynamic modelaii^. The presence of two types of 
charges leads to the situation that the boundary condi- 
tions of Maxwell's equations are not sufficient to solve 
the problem. The problem has been originally identified 
by Pekar— , with a variety of "additional boundary con- 
ditions" (ABC) proposed over the years^^Ti^i, and com- 
pared with experimental data3. In Ref. |32| . a model with 
two charge types, similar to ours, was considered and an 
ABC was proposed by arguing that dissipated energy is 
conserved across a boundary. It is known, however, from 
ABC studies22r— that microscopic information must be 
used, i.e. Maxwell's equations by themselves are insuffi- 
cient. 

We address here the ABC problem by solving a Boltz- 
mann equation, allowing for impurity scattering in the 
bulk as well as surface scattering that mixes CC and CL. 
Note also that the CL allows charge conservation to hold 
even if the bulk current perpendicular to the surface is 
locally finite, representing deviations from specular re- 
flection. In this sense the localized charge is a measure of 
diffuse scattering or disorder at the surface: the scattered 
bulk current becomes a surface current which is allowed 
to diffuse along the surface. A key point of our model are 
the transition rates between the two types of charges that 
provide a restoring force (plasma oscillation) and damp- 
ing for the combined charge oscillations. By allowing 
for bulk and surface impurity scattering, our approach 
complements traditional band structure theoriesi^i^ and 
provides a unifled framework for analyzing SP dispersion 
and linewidth. The resulting requations involve a few pa- 
rameters of the poorly known CL and of the surface-bulk 
charge scattering. We start with discussing the qualita- 
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tive effect of the surface scattering and then derive the 
actual SP dispersion, showing a variety of forms corre- 
sponding to experimental data. 

Qualitative features. Our model implies that the 
surface scattering that mixes CC and CL enhances the SP 
linewidth F. Hence F is relatively small if the charge re- 
sponse is dominantly CC (large positive slope, B > Cki 
with ki = 1 A being a typical large momentum) or 
dominantly CL (large negative slope, B < — Cfci, or 
small slope, \B\ < O.lC/ci). Table U shows aU the corre- 
sponding cases that we are aware of, confirming this trend 
for small damping. Cases where O.lC/ci < \B\ < Cki 
may or may not correspond to dominant CC or CL and 
a detailed fit oiuj{k) and F(fc) is needed to determine the 
parameters. 

The model and nonlocal electrodynamics. We 

proceed to define the charges and currents of the model. 
The bulk medium is located in z < 0, and the surface 
layer in < z < d. We write p{r) and Ps(r), respectively, 
for the charge density in the two regions. We expand the 
current response, considering small charge gradients 



-j + dtj 

T 

is + dtjs 



in 



E - clVps 



z<0 



< z < d. 
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The equilibrium charge density provides restoring forces 
LOp/An and lu'^/Att in the bulk and surface regions, where 
ujp^s are the corresponding plasma frequencies if the 
CC/CL were decoupled. The bulk and surface compress- 
ibilities are expressed in terms of the sound velocities 
Cb^s, respectively, while c^r and c^Tj are the correspond- 
ing diffusion coefficients. An important realization of a 
surface layer is the case of a charge spill-out, as in the 
alkali metals^iiiiS. In this case, the CL is further away 
from the restoring force provided by the equilibrium situ- 
ation. Hence the equilibrium charge, when averaged over 
the CL thickness d in the oscillating state, is reduced and 
the restoring force is weak, ujs < i-^p- Note that even with 
Ws = 0, there is eventually a restoring force via coupling 
to the bulk so that cj(fc — >■ 0) = ujp/V^. We also note 
that when LUg — 0, Feibelman's sum rule Eq.([T|)^^ is still 
valid with z = as origin, since the charges at z > 
respond only in a nonlocal way via the gradient Vp^ 
which affects higher order 0{k^) terms. When Wg ^ 0, 
e.g. for nanostructures with their own ions and equilib- 
rium charge density, the sum rule ([U is modified such 
that 



B = 



2^/2 



J dzz5n{z) 
J dz Sn(z) uj, 
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adding a positive term to the slope. 

In the following we omit the damping terms (1/t, 1 /ts) 
for brevity. They are easily restored by multiplying the 
bulk [surface] parameters ojp and [oj^ and c^] with the 
factor 1 -I- i/{uJT) [1 + i/{ujTs)], respectively. 



We look for solutions to the fields at frequency ui 
that vary with the wave vector k parallel to the surface. 
Charge conservation then determines the charge profiles 



z < : 
< z < d : 



4^' 



-Lu^]p{r) = 



(4) 



p(A:,z) = po(fc)e'''^ 
[c.2-c2v2-c.2]p,(r) = 
Ps(fc, z) — ja{k) coshvsZ + 7i(fc) sinhws-z 



The inverse decay lengths Vb and Vs for the bulk and 
surface charges are 



Vs = — + c2fc2 - LU^ , 



Vb = 



Cb 



oj?, + clk"^ — 



(5) 



Current conservation at z = d yields {ljI / A:n)Ez{k,d) 
c^dzPsik, d) = 0, while at z = 0, 



^i?,(fc,0- 



cld:,p{k,0~ 



= -^Ez{k,0+) - cldzps{k,0+) (6) 

where 0^ denotes the limit z — )■ ±0. 

Maxwell's equation determines the longitudinal com- 
ponents of E in terms of the charges, however the trans- 
verse components require four unknowns: two ampli- 
tudes for upward and downward propagation in the layer, 
and one amplitude each in bulk (z < 0) and in vacuum 
(z > d). Together with the amplitudes /3o,70j7i for the 
charge density, we have seven unknowns. Each interface 
provides three relations (the matching of Ez, dzEz and 
jz), hence one boundary condition is missing. The ne- 
cessity for an additional boundary condition (ABC) has 
a long history and occurs when more than one mate- 
rial mode is present, for example the 70, 71 modes of Eq. 
(|4p2r— . Therefore a microscopic input is needed to pro- 
vide the ABC. 

Before developing the ABC, we discuss some general 
properties of the CC-I-CL system. At very small parallel 
momentum, k < cup/c ~ 2 • 10~^A , the SP disper- 
sion approaches the light dispersion uj{k) -> ck. While 
this range can be optically probed^'^, most of the data 
is found by electron energy loss spectroscopy at higher 
k values. Therefore, we focus on the non-retarded re- 
gion uip/c <^ k, formally taking c — >■ cx). Since the limit 
uj{k — >■ 0) = ujp/ is a bulk property, we expect that the 
parameter Ws has a small effect. In fact, the numerical 
solutions shown below confirm this for ujs ^ 0.3 Up and 
k < OAujp/cb- Hence we display the result for ojg — 0, 
which has a relatively simple form: 



^Po(fc)[w6 



Lo'^Vb 



-kd 
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coshwsd 



l)-i;,(2— -l)tanh(i;,d)]. 
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Ref. 


sample 


A (eV) 


B (eV-A) 


C (eV-A ) 


r (eV) 


6 


K 


2.73 


-0.95 


0.09 


0.3 


4 


Ag(llO) k II [001] 


3.76 


- 


3.68 


0.11 


4 


Ag(llO) k II [110] 


3.86 


- 


1.46 


0.18 


4 


Ag(lll) 


3.69 


- 


4.17 


0.06 


9 


0.05ML O/Ag(001) 


3.71 


- 


3.1 


~ 0.3 


9 


O.IML O/Ag(001) 


3.70 


1.1 


-0.9 


~ 0.3 


9 


0.15ML O2/Ag(001) 


3.69 


0.86 


-0.1 


~ 0.3 


11 


Ag(OOl) 


3.71 


1.45 


- 


0.15 


13 


Cu(lll) 


1.18 


-0.72 


- 


0.2 


14 


lOML Ag/Ni(lll) 


3.75 




1.57 


~ 0.1 


16 


sputtered 22ML Ag/Cu(lll) 


3.76 


-0.08 


2.50 


0.18 


7 


Mg(OOOl) 


7.38 


-3.023 


9.78 


1.2 


10 


Al(lll) 


10.86 


-3.1 


7.7 


2.3 



TABLE L Data on SP dispersion tj = A + Bk + Cfc^ and linewidth F at small k. The coefficients yl,_B,C are given by the 
corresponding references, with "-" indicating a coefficient that was assumed to be zero for the fit. All cases with |_Bj > Ck\ 
or |_B| < O.lCfci are shown with a relatively small F (momentum scale fci = lA ^). The last two entries show cases where 
O.lCfci < 1B| < Cfci with a large F. 



where the ABC is needed to fix the ratio po (fc) /7o (fc) • For 
a CC system (70 = 0) we recover the well known^''^'^^ SP 
dispersion = \uj^ + iS^kjlvh with a dominant linear 
dispersion at small /c, i.e. 

CC: c.(fc) = ^ + £^+0(fc2) (8) 

For the CL system (po = 0), we obtain to lowest order 
as d — >■ a purely quadratic dispersion uP' = + c^fc^, 
featuring the surface speed of sound Cg. Keeping d finite, 
and expanding for small fc, we get 

CL: w(fc) = i::J^-£i^tan(d72) + 0(fc2), (9) 
V 2 2 

where d — ujpd/ {■\/2cs)- Note the negative slope, con- 
sistent with the charge being outside the bulk ions (at 
z > O)'^^'. The divergence at tan J/ 2 = 00 corre- 
sponds to a resonant multimode response^ii^i^, which 
may appear at higher frequencies if the thickness d is 
large enough. Indeed, for cj^ = as assumed here, 
we find that this corresponds to a vanishing total layer 
charge, dz ps{k, z) = 0: these are surface dipoles or 
higher multipoles. They have been observed in alkali 
metals^'^'^° where the charge response outside the bulk 
is significant. These modes, having larger linewidths, are 
not considered in the solutions below. 

We note also that in order to get an acoustic plas- 
mon solution^'^^ with Ld{k) — > as fc — )■ 0, one would 
require in Eq.(I7]) jod — > —po/vh corresponding to oppo- 
site signs of the total bulk {po/vt) and total surface (70^) 
charges. A specific ABC is then needed to achieve acous- 
tic plasmons-. Our model allows for surface scattering 
between the two types of charges, a distinct situation. 



Additional boundary condition via Boltzmann 
equation. We proceed to derive the ABC by us- 
ing a Boltzmann transport equation for the charge 
response^^^. This approach is known to be notori- 
ously difficult when matching conditions at an inter- 
face are involved^. We have adopted an alternative ap- 
proach and consider a bulk half-space terminated by a 
two-dimensional surface sheet that represents the charge 
integrated over the layer < z < d. The matching 
conditions are then replaced by surface scattering that 
mixes bulk and surface charges. At equilibrium, the bulk 
electron distribution per phase space d^xdzd^p/ (27r/i)"^ is 
/"^(x, z, p) while the surface distribution per phase space 
d^xd^q/ {27Th)^ is /°(x, q). The surface charge sheet in- 
volves only the lateral coordinates x = {x,y) and a two- 
dimensional momentum q. We note that this phase space 
description is possible even if the surface states are local- 
ized, by working with the Wigner transform'^^''^^. In re- 
sponse to a weak electric field, the distributions become 
/° + /, f° + fs, with /, /s < 1. The transition probability 
due to bulk impurity scattering per d'^pdt from state p 
to p' is Vl^(p',p). Within the Born approximation, or 
more generally in the presence of time reversal and space 
inversion, we have M^(p,p') = W{p',p). To 1st order in 
/, E, with a bulk velocity v(p) 

[dt + v(p) • V]/(x, z, p) + eE • Vp/" (x, z, p) 

= / Ty(p',p)[/(x,z,p')-/(x,z,p)] (10) 
Jp' 

with the shorthand J^, = Jd^p' /{2TTh)^ . At the sur- 
face we have a scattering cross section W^' (q', q) between 
surface states, as well as scattering of bulk states with 
Pz > into the surface {Ws{(l, p)) and scattering of sur- 
face states back into bulk states with p^ < {Ws{p, q)). 
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As in the bulk, we assume W{q',q) = M^'(q, q'). For 
W^s(ctj p) we can assume reflection symmetry only in x, y 
as well as time reversal, hence Ws{c[, p) = Ws{—c[, — p) = 
Ws{p,<l), where p = {px,Py,—pz)- Boltzmann's equa- 
tion, to first order in E, /, fs becomes, with a surface 
velocity u(q) 

[dt + u(q) • V]/,(x, q) + eE • Vq/,"(x, q) 

= / W^'(q',q)[/.(x,q')-/.(x,q)] 

+ f W^,(q,p){/(x,0,p)-/,(x,q) 

+ /0(x,q)[/(x,0,p)-/(x,0,p)]} (11) 

Note that unlike the bulk equation (ITO)) . /" does not can- 
cel. The final ingredient is a matching condition: the to- 
tal flux with velocity ±Vz is matched with the probability 
of scattering in or out of surface states. Taking Vz > 
and considering the Ws(q, p) symmetries, 

t>,[/(x,0,p)-/(x,0,p)] = 

/ t^.(q,p){/(x,0,p)-A(x,q) 

+ /^(x, q)[/(x, 0, p) - /(x, 0, p)]} (12) 

This relation assumes specular reflection in the absence 
of surface states, while the latter produce nonspecular 
reflection in each channel. Integrating Eq. p^ over 
d'^p, Pz > yields the total current incoming into the sur- 
face, j™(x), and by comparing with Eg. pTj) integrated 
over d^g, we get charge conservation at the surface 

aiP,(x) + v.j.(x)=jT(x) (13) 

where p^(x) = J^fs{x,q), js(x) = u(q)/,(x, q). 

To solve the transport problem, we make the Pi ap- 
proximation, i.e., on spatial scales large compared to the 
mean free path, only the weakest angular dependence in 
the momentum variable is maintained^. We therefore 
use the following expansion 



/(x,2;,p) 
/s(x,q) 



p(x, z)/(i)(p) + ^ j,(x, z)v,{p)gll\p) 



p,(x)/ii)(g)+j,(x).u(q)5W(g) (14) 



where p, q denote the modulus of the vectors p, q. The 
reduced symmetry of the surface allows for a diagonal 
tensor with elements gi]^{p) = gyy{p) 7^ g^zz{p)- For 
elastic scattering, we expect the dependence on p, q to 
be localized near the corresponding Fermi surface, if the 
latter is well-deflned (the surface states may be spatially 
localized). 

Integrating the v or u moment of Eqs. (jlOllip 
yields Eq.Q, identifying the plasma frequency as 
— Ane^ ^/''(x, z, p) (assuming diagonal response) 

which becomes ojp = Aime^ /m* for a parabolic dispersion 
with an effective mass m* . Similarly = Annsc'^ /m*d 



where n, are equilibrium bulk and surface densities, 
respectively. The v or u moments of Eas. (|10llip yield 
the sound velocity — vf f^^\p) and similarly for Cg, 
and also identify the bulk and surface relaxation times 



W^(p',p)«f(p)5l''(p) 



/ I^'(q',q)«,2(q)g(i)(q) 
J a. a' 

W.M,p)u^,{ci)g^\q) (15) 



1 

T 
1 



p,Pi>0 

where 1/r may in principle be anisotropic. We also note 
that the moment of Ea.(|12p requires in the expan- 
sion (|14|) of the momentum distribution a quadrupole 
term /(x, z, p) ~ v^Vz- This term contributes to shear 
viscosity at the surface. It is of higher order in veloci- 
ties, however, and does not affect the Vz moments that 
we discuss now. 

Taking the Vz moment of the matching condition ()12|) , 
we encounter the following integrals: 



Ws (q, p)vz [p(x, 0)/(i) {p) + Vzjz (x, 0)g W (p)] 



p,p^>0 J q 

a2p(x,0) -t- a3j2(x,0) 



/ / Ws{q,p)vzfi'Hq) = ai 
/ / W^.(q,p)2^2^W(p)/r(x,q) 

p,Ps>0 q 



"5 



(16) 



Note that the difference — as involves an integral over 
1 — 2/"(x, q), which vanishes if the surface states have 
an electron-hole symmetry. 

Putting everything together, the Vz moment of Eg. p^ 
becomes the ABC 



r(x) 


= ac6po(x) 




a2 


acb 






ai -1- q;3 ■ 


1 


04 


7 


ai-\- 



"5 



"5 



(17) 

(18) 
(19) 



The two coefficients on the rhs can be interpreted as a 
probability of a bulk carrier getting trapped in a surface 
state (coefficient a) and a desorption rate back into the 
bulk (rate 1/t'). 

Solutions for the dispersion relation. To make 
contact with the electromagnetic calculation, we iden- 
tify ps with the charge density averaged over the layer 
< z < d, i.e., integrating Eq.®, we get Ps{k) — 
[7o(fc) sinhugrf — 7i(^)(l — coshvsd)]/vs. When the 
boundary condition for the electric field Ez and Eq. ([6]) for 
the current flowing into the layer, are combined with 
the ABC (fTTl) . one gets a linear system for the charge 
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0.1 0.2 0.3 0.4 cOp 0.1 0.2 0.3 0.4 ojp 



FIG. 1: {left) Surface plasmon dispersion relation uj{k) for different values of the surface scattering parameter a in Ea. p7|l : 
a = 0.1 (thin blue lines), ol — 1 (dashed red lines) and a = 10 (thick green lines). Upper lines starting at a;p/%/2: real part; 
lower lines: imaginary part x(— 1). Other parameters: uja = 0, Ca/cb = 1, ojpT = 10, ujpTs — 10, t' — oo and d = ujpd/cb ~ 3.5. 
(right) Effects of a nonzero restoring force in the surface layer (plasma frequency ujs): ujs — (thin blue curves), tOa/i^p = 0.5 
(dashed red curves) and LJa/oJp = 0.9 (thick green curves). The upper curves are the real part and the lower ones are the 
imaginary part x(— 5). The other parameters are Ca/cb = y/S, oOpT = 10, oOpTa = 10, r'up = 1, a = 0.1 and d = 7. 



amplitudes pQ{k),jo{k),ji{k). For the interesting case 
LUs — 0, this becomes 



acbPoik) = [i 



1 , tanhiigd 



70 (fc) (20) 



and with Eq.(l7l) yields the SP dispersion 



Cb- 



U6(2w2 



ulk 



Vsk{ 



1 



sinhucd tanhw,(i 



)-vl{2- 



(21) 



The parameter a thus measures the relative importance 
of the CC and CL amplitudes. For a ^ 1, the dispersion 
corresponds to a CC oscillation, Eq.(|Sl), vifhile for a 1 
the dispersion is dominated by the CL, Eq.([5]). We recall 
that a large value of a is facilitated when electron-hole 
symmetry holds [a^ = as in Eg. dTSl) ]. 

We display in Fig IT] numerical solutions for the SP dis- 
persion. On the left the effect of the CC/CL mixing 
parameter a is shown: the linewidth is relatively small 
for a — 0.1 (dominant CC) or a = 10 (dominant CL), 
while it is much larger for a — \ corresponding to strong 
CC/CL mixing. This figure is for d = uj.pd/cb — 3.5 
demonstrating also a negative initial slope when the CL 
is significant, i.e. a = 1, 10. When d is smaller, main- 
taining a > 1, the linear term is suppressed (Eq. [T|) and 
the dispersion becomes purely parabolic, as observed in 
some cases^ii^. Fig[TJright) demonstrates the effect of 
ojs and the possibility of generating a local minimum in 
r(fc), as observed in some cases^^-'- ^. This situation is 
facilitated by uJs ~ ujp and r'ujp = l,d = 7 . The effect of 
ujs on the real part is to increase the linear term in k, as 
discussed in Eq.(|31). 



Discussion. We developed here a general scheme for 
solving the SP dispersion at low fc, allowing for two 
charge types, CC and CL, both with nonlocal diffusive 
response. We have solved the notorious ABC problem 
by a Boltzmann equation that allows for CC/CL sur- 
face scattering. The resulting dispersion is consistent 
with a large variety of data accumulated on pure met- 
als (where the CL corresponds to the spill-out charge), 
on sputtered systems, disordered surfaces, quantum well 
surface states, all leading to distinct electron states in a 
surface layer. 

We discuss now the role of the most significant pa- 
rameters: surface mixing a, layer thickness d, surface 
restoring force tol , and how they correspond to observed 
data. Consider first the surface mixing: when it is weak 
(small or large a) , it leads to a relatively small linewidth 
as in the first eleven lines of Table HI Weak surface mix- 
ing corresponds also to dominant CC or CL which is re- 
alized e.g. when the dispersion is dominantly linear or 
quadratic, as indeed in those eleven lines. Cases where 
the linear and quadratic terms are comparable may (last 
two lines of Table H]) or may not lead to a large linewidth. 
In the latter case a detailed fit is necessary, and is also 
demonstrated in Fig. [TJleft). 

The layer thickness d allows for a negative linear dis- 
persion, as seen in many metals (Table |l|. When d — > 
and CL dominates (a ^ 1), the dispersion is purely 
quadratic, as indeed seen in several cases^i^ii^. 

The surface restoring force Ws is expected to be small, 
LUs < ujp, in the case of alkali metals where the spill- 
out charge is further away from the equilibrium charges. 
In other large lUs shows the peculiar feature 

of a minimum in r(fc), as demonstrated in Figlljright). 
This minimum, as Al(lll)-1" and Ag on Si(lll)^^, has 
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been discussed in terms of specific band structures. Our 
approach is an alternative description in terms of the 
charge distribution. With a few parameters our model is 
thus able to describe SP dispersion and linewidth on a 
large variety of metals and surfaces, as seen over the last 
decades. 

Acknowledgments. We are thankful for stimulating 
discussions with E. V. Chulkov, D. Davidov, O. Entin- 



Wohlman, M. Golosovsky, F. Guinea, and V. M. Silkin. 
We also thank E. Eizner for significant help with data 
analysis. This research was supported by a Grant from 
the G.I.F., the German-Israeli Foundation for Scientific 
Research and Development. 



^ F. Garci'a-Moliner and F. Flores, Introduction to the theory 
of solid surfaces, Cambridge University Press (1979). 

^ For a review see: J. M. Pitarke, V. M. Silkin, E. V. 
Chulkov, and P. M. Echenique, Rep. Prog. Phys. 70, 1 
(2007). 

^ J. Homola, Ed., Surface plasmon resonance based sensors, 
Springer (2006). 

" S. Sutto, K.-D. Tsuei, E. E. Plummer and E. Burstein, 
Phys. Rev. Lett. 63, 2590 (1989). 

^ M. Rocca and V. Valbusa, Phys. Rev. Lett. 64, 2398 (1990) 

® K. -D. Tsuei, E. W. Plummer, A. Liebsch, E. Pehlke, K. 
Kempa and P. Bakshi, Surf. Sci. 247, 302 (1991). 

^ P. T. Sprunger, G. M. Watson and E. W. Plummer, Surf. 
Sci. 269-270, 551 (1992). 

® M. Rocca, Li Yibing, F. Buaticr dc Mongcot and U. Val- 
busa, Phys. Rev. B 52, 14 947 (1995). 

® L. Savio, L. Vattuone and M. Rocca, Phys. Rev. B 61, 
7325 (2000) 

G. Chiarello, V. Formoso, A. SantanicUo. E. Colavita and 

L. Papagno, Phys. Rev. B 62, 12 676 (2000). 

L. Savio, L. Vattuone and M. Rocca, Phys. Rev. B 67, 

045406 (2003). 

12 Y. Yu, Y. Jiang, Z. Tang, Q. Guo, J. Jia, Q. Xue, K. Wu 

and E. Wang, Phys. Rev. B 72, 205405 (2005). 
1^ A. Politano, G. Chiarello, V. Formoso, R. G. Agostino and 

E. Colavita, Phys. Rev. B 74, 081401 (2006). 
A. Politano, R. G. Agostino, E. Colavita, V. Formoso and 
G. Chiarello, phys. stat. sol. (RRL) 2, 86 (2008) 
1^ A. Politano, G. Chiarello, V. Formoso, Plasmonics 3, 165 
(2008) 

A. Politano and G. Chiarello, J. Phys. D 43, 085302 (2010). 
1^ A. Varykhalov, A. M. Shikin, W. Gudat, P. Moras, C. 
Grazioli, C. Carbone, and O. Rader, Phys. Rev. Lett. 95, 

247601 (2005). 

1® A. Liebsch, Phys. Rev. Lett. 71, 145 (1993); Phys. Rev. B 

15, 11317 (1993). 
1® B. Diaconescu, K. Pohl, L. Vattuone, L. Savio, P. Hof- 

mann, V. M. Silkin, J. M. Pitarke, E. V. Chulkov, P. 

M. Echenique, D. Farias, and M. Rocca, Nature 448, 57 

(2007). 



2° Q. A. Turchette, D. Kielpinski, B. E. King, D. Leibfried, 

D. M. Mcckhof, C. J. Myatt, M. A. Rowo, C. A. Sackctt, 

C. S. Wood, W. M. Itano, C. Monroe, and D. J. Wineland, 
Phys. Rev. A 61, 063418 (2000). 

D. Leibrandt, B. Yurke, and R. Slusher, Quantum Inf. 
Comput. 7, 52 (2007). 

22 C. Henkel and B. Horovitz, Phys. Rev. A 78, 042902 

(2008). 

2^ R. Dubcssy, T. Coudreau, and L. Guidoni, Phys. Rev. A 

80, 031402(R) (2009). 
2* N. Daniilidis, S. Narayanan, S. A. MoUer, R. Clark, T. E. 

Lee, P. J. Leek, A. Wallraff, S. Schulz, F. Schmidt-Kaler, 

and H. Haffner, New J. Phys. 13, 013032 (2011). 
25 R. H. Ritchie, Surf. Sci. 34, 1 (1973). 
2« P. J. Fcibclman, Phys. Rev. B 40, 2752 (1989). 
2^ S. I. Pckar, Zh. Eksp. Tcor. Fiz. 33, 1022 (1957) [Soviet 

Phys. JETP 6, 785 (1958)]. 
2* R. Zeyher, J. L. Birman and W. Brenig, Phys. Rev. B 6, 

4613 (1972). 

2^ C. Schwartz and W. L. Schaich, Phys. Rev. B 26, 7008 
(1982). 

3° H. C. Schneider, F. Jahnke, S. W. Koch, J. Tignon, T. 

Hasche, and D. S. Chemla, Phys. Rev. B 63, 042202 (2001). 

For a recent critique and references see M. G. Silveirinha, 

New J. Phys. 11, 113016 (2009). 
''*2 F. Forstmann and H. Stenschkc, Phys. Rev. Lett. 38, 1365 

(1977); Phys. Rev. B 17, 1489 (1978). 

V. Lirtsman, R. Ziblat, M. Golosovsky, D. Davidov, R. 

Pogreb, V. Sacks-Granek, and J. Rishpon, J. Appl. Phys. 

98, 093506 (2005). 
3* A. J. Bennet, Phys. Rev. B 1, 203 (1970). 
^5 N. W. Ashcroft and N. D. Mermin, Solid State Physics, 

Harcourt College Publishers (1976), Ch. 16. 

J. M. Ziman, Electrons and Phonons, Oxford University 

Press (1960), Ch. VII, XL 
^"^ T. R. Kirkpatrick, Phys. Rev. B 33, 780 (1986). 

C. Henkel, C. R. Acad. Sci. (Paris), Serie IV 2, 573 (2001). 
3^ G. G. M. Coppa and A. D'Angola, Eur. Phys. J. D 6, 533 

(1999). 



